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We propose that all light fermionic degrees of freedom, including the Standard Model (SM) 
fermions and all possible light beyond-the-standard-model fields, are chiral with respect to some 
spontaneously broken abelian gauge symmetry. Hypercharge, for example, plays this role for the 
SM fermions. We introduce a new symmetry, U{l)v, for all new light fermionic states. Anomaly 
cancellations mandate the existence of several new fermion fields with nontrivial U{l)i, charges. We 
develop a concrete model of this type, for which we show that (i) some fermions remain massless 
after U{l)i, breaking - similar to SM neutrinos - and (ii) accidental global symmetries translate into 
stable massive particles - similar to SM protons. These ingredients provide a solution to the dark 
matter and neutrino mass puzzles assuming one also postulates the existence of heavy degrees of 
freedom that act as “mediators” between the two sectors. The neutrino mass mechanism described 
here leads to parametrically small Dirac neutrino masses, and the model also requires the existence 
of at least four Dirac sterile neutrinos. Finally, we describe a general technique to write down 
chiral-fermions-only models that are at least anomaly-free under a U(l) gauge symmetry. 


I. INTRODUCTION 

All confirmed fundamental fermion fields are chiral. We can express all known matter fields as left-handed Weyl 
fermions - Q, d‘^, L, e’^ - that transform under the Standard Model (SM) gauge group SU{3) x SU{2) x U(l) 
as ( 3 , 2 )i/g, ( 3 , l)_2/3, ( 3 , 1)1/3, ( 1 ) ‘^)-i /2 and ( 1 , l)i, respectively. There are no vector-like fermions, i.e., pairs 
of fields with ‘equal-but-opposite’ charges. One consequence of this experimental fact is that, before electroweak 
symmetry breaking (EWSB), all fundamental fermions are massless. A corollary is that all fermion masses are 
proportional to the parameter that controls EWSB. In the case of the SM, this is the vacuum expectation value v of 
the neutral component of the scalar Higgs field H, ( 1 , 2)1/2. 

One can speculate why, to date, all identified fermion fields turned out to be chiral. On one hand, the existence 
of a single chiral fermion charged under a gauge symmetry implies the existence of several other chiral fields in order 
to cancel the gauge anomalies. On the other hand, vector-like fermions, if they were to exist, would have masses 
unrelated to EWSB and could be out of the reach of current experimental probes. The lack of concrete, unambiguous, 
experimental evidence leads one to conclude that, if they do exist, vector-like fermions are either very heavy or, in 
the event that the new fermion is not charged under the SM gauge group, very weakly coupled. 

Yet, we know the SM is not complete. The SM fermion content, combined with the SM Higgs sector, leads to the 
prediction that neutrino masses are zero. While this is not a bad approximation -- neutrino masses are known to be 
tiny - it is factually incorrect M. The origin of nonzero neutrino masses is currently unknown, but it is clear that 
new degrees of freedom must exist. The preferred neutrino mass mechanism makes use of the fact that neutrinos are 
singlets of the unbroken SMgauge groups and hence can have Majorana masses. In the simplest realizations, the so 
called Seesaw Mechanisms (a-Qj ^ heavy state is introduced that mediates the Weinberg operator, {LH){LH). On 
the other hand, Dirac neutrino masses also require the introduction of a new, SM gauge-singlet fermionic degree of 
freedom. While technically natural, this hypothesis is often considered to be contrived since it requires extremely 
small Yukawa couplings. 

New degrees of freedom are also required in order to explain the Dark Matter (DM) puzzle. Regardless of its 
nature, the existence of DM also appears to imply a new fundamental physics scale (see Q for an attempt at a 
counter example). 

In this paper, we propose that all light fermionic degrees of freedom are chiral in the sense described above - 
what we mean by “light” will be made clear below. In particular, new light fermionic degrees of freedom required to 
explain neutrino masses and DM must also be chiral. To achieve this, we postulate the existence of, at least, a new 
nonanomalous C/(l) gauge symmetry that we dub 17(1)^. We impose, similar to the SM, that (i) all new fermions are 
charged under U{l)i, - that plays the role of SM hypercharge - and (ii) no light vector-like fermions exist. Anomaly 
cancellations then require the existence of several new fields. With these ingredients we construct a model for a Dark 
Sector (DS) that, in analogy to the SM, shows accidental global symmetries as a consequence of the charge 

assignments. These in turn imply the existence of at least one stable massive particle that serves as a DM candidate. 

Similar to the SM active neutrino, massless fermions appear in the DS. We show that the smallness of neutrino 
masses can be understood if the initially massless DS states are charged under the SM lepton number symmetry. 
Lepton number is communicated between the SM and the DS via a nonchiral, heavy mediator sector charged under 
both the SM and the DS gauge symmetries. After integrating out the mediator sector, this model produces naturally 
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small Dirac neutrino masses through a mechanism first described in Q. This contrasts with previous models of a 
chiral DS, e.g. , where the Seesaw is still invoked. 

Kinetic mixing between the SM and the gauge bosons [llj is unavoidable in this scenario. It leads to several 

consequences, most importantly, the possibility of detecting the DM particles in the laboratory [I^. Barring the 
possibility of fine tuning, the predicted DM cross section would be within reach for the next generation of DM direct 
detection experiments. 

In order to develop the scenario described above, an algorithm is necessary to write chiral models that are at least 
anomaly-free under a C/(l) gauge symmetry [I^ - see also |l4fl5l|. A general technique to do so is described in Sec. lVIi 
towards the end of this manuscript. Indeed, following the steps detailed in Sec. IVIi one can construct any number of 
models that fit our requirements. The method also allows one to address several technical questions including what 
are the “minimal” anomaly-free U{1) gauge theories with chiral fermions. This question is answered precisely, for two 
different criteria for minimality: smallest highest charge and smallest number of fermion fields. 

Before embarking in such general considerations, we describe in Section|TT]a simple yet phenomenologically appealing 
realization of our scenario and, in Section lllll discuss in some detail some of the relevant features of this model, mostly 
those related to the introduction of new gauge interactions. We discuss how small neutrino masses and dark matter 
can be accommodated in this model in Sections m and |Vl Finally, in Section IVIIi we briefly discuss other models, 
potential research directions, and provide some concluding thoughts. 


II. CHIRAL U{1)^ MODEL 

In this section we construct a Lagrangian consisting of the SM, a I/(l),y-charged, nonanomalous “dark sector (DS)” 
and a “mediator sector”. 


— -SfsM + -Sfbs + -SfMix + -^Med, (1) 

where and ^ued represent the SM, DS, and Mediator Lagrangians, respectively. The term contains 

renormalizable operators that “mix” the SM and DS degrees of freedom, including the kinetic mixing of with 

the hypercharge U{1)y and the scalar potential coupling between the Higgs field and the equivalent “dark Higgs” 
scalar field. We discuss these term in detail in Sec. nni 

In addition to the SM gauge group, .if is also invariant under a gauged U{l)i, symmetry, which “lives” in the DS. 
We assume the fermionic particle content, assuming all fermions to be chiral, is as follows, keeping in mind that all 
fermions are left-handed Weyl fields: 

• three fields with charge -fl - 1+^’^; 

• two fields with charge —2 - 2^’^; 

• two fields with charge —3 - 3^’^; 

• three fields with charge -1-4-4^^’^; 

• one field with charge —5 - 

We adopt a normalization for the new gauge coupling g^, in which all t/(l)j/ charges are integers. The SM fields are 
not charged under [/(l)iy, while the new fermions defined above are not charged under the SM gauge symmetry.* 
f7(l)i/ is assumed to be spontaneously broken in order to render most of the new fermions and the new gauge boson 
massive. We achieve this by, similar to the SM, adding one scalar field (j) with charge -1-1 and a scalar potential such 
that (j) has a nonzero vacuum expectation value v^. 

Having defined the gauge symmetry and the particle content, the renormalizable Dark Sector Lagrangian is well- 
defined: 


-Sfos = -Sfos-kin + ^DS-Yuk + 


( 2 ) 


where .ifos-kin and -Sfos-Yuk represent the kinetic-energy and Yukawa terms respectively, while V{(j)) is the scalar 
potential for the field </>. The kinetic-energy terms are 
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^DS-kin = 




( 3 ) 


From the point of view of the SM, all new fermions are gauge-singlet “neutrinos,” hence the name U(l)iy 
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where 


, ( 4 ) 

is the U{l)u field strength and B^ the U{l)i, gauge field. j = 0, 1, 2, fc = l, 2. Note that, here, 1+ (which has 
charge —1) represents the complex-conjugated, right-handed Weyl field associate to the 1+ left-handed Weyl field, 
etc. 

The Yukawa interactions in the DS are given by: 

- i?DS-Yuk = + hM + h,k + h.c., (5) 

where /, h are 2x3 and 3x3 Yukawa matrices respectively. In what follows the indices i, k will be omitted whenever 
there can be no confusion. 

After U{l)i, breaking, the 3_, 5_ and 4+ fields “pair up” into three Dirac fermions, labelled from here on y*, their 
masses-squared associated to the eigenvalues of the 3x3 mass-squared matrix where 

= h^jV^ , i, j = 0, 1, 2. (6) 

In general, is expected to have three nonvanishing eigenvalues and hence describe three massive Dirac fermions. 

Similarly, the two 2_ fields pair up with two linear combinations of the three 1+ fields into two Dirac fermions that 
we label N^. Note that there is no 2^ field. Hence, there is a massless chiral linear combination of 1+ fields, that 
we call 2 ^“. All the masses-squared, including the vanishing one, can be obtained from the eigenvalues of the 3x3 
mass-squared matrix where 


I b J — 0, 1, 2 , (7) 

and where we define /^o = 0 for all * = 0,1, 2. In summary, in the absence of fields that connect the DS to the 
SM, after symmetry breaking, the model contains 3 4-2 = 5 massive Dirac fermions and one massless 

left-handed Weyl fermion, 

The unitary rotation of the fields 1+ and 2_ that renders diagonal has no effect in Eq. (jS]) since these rotations 
operate independently on fields of the same charge. Performing a similar change of basis on the 3_, 4+ and 5„ fields 
in order to render diagonal is different since, in this case, the 5^ and 3^;^ fields must be grouped together as 
components of a 3-vector. If is diagonalized by the transformation 

UM^V^ = diag{M^,, } , ( 8 ) 

where U, V are 3x3 unitary matrices, “mixing” appears among the x fields. The U (1)^ couplings to the Xi mass eigen¬ 
states are off-diagonal which indicates the presence of “flavor-changing neutral currents” within each DS generation. 
Explicitly, in the mass basis, the B couplings are given byi 

gu - NrXiiNr - 2NLXf^NL - ixRlfj-XR + Xl2357mXl) (9) 

where the L, R subindices mean that we take the left- or right-handed components of the corresponding Dirac 
field. Indices running over fields of the same charge are implied and the matrix Q 35 in the last term is the rotated 
charge-matrix 


235 = V 




( 10 ) 


Since this matrix is nondiagonal, generically, two of the x fermions in the 345-sector, say X 2 and xa, are unstable; 
decaying for example into a lighter x field and a pair via B boson exchange. This B boson can be on- or 

off-shell depending on the masses of the Xi- Also notice that the stable, lightest y-particle, xi, couples to B with a 
right-handed coupling —Ag^ and a left-handed coupling (—3 -I- (—5 -I- 3 )|V^ 5 p) 5 iy, where |V^ 5 p is the probability that 


t The Dirac fermions are (ignoring i,k indexes), N'^ = (2_, 1+), = (3—,4+), (5_,4+), where the comma indicates that x are linear 

combinations of those two Dirac fermions. 
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a xi state will interact as a charge —5 object. The scattering cross-section of an unpolarized xi beam via B exchange 
will be proportional to an effective charge-squared 


Qx = ^(4'+[(3 + 2|Vx5|"]')e[12.5,17], 


( 11 ) 


which will come in handy later. 

After all the smoke has cleared, Jfos contains four accidental global 17(1) symmetries, 17(1)x 17(1)jvi x 17(1)jv 2 x 
17(l)i/c and hence four stable particles. More structure in the DS will, in general, reduce the accidental global 
symmetry of the Lagrangian. For example, if there is a second 17(1) i, charged scalar field (j)', also with charge -1-1, the 
17 ( 1 )'^ global symmetry would be, in general, reduced to 17 ( 1 )-^ x 17(l)i2, where 12 refers to the fields with charges 
-1-1, —2. In this case, only the lightest x held would be absolutely stable, along with the massless v'^ state. Ultimately, 
we will associate the lightest x field with dark matter, while will play the role of the left-handed antineutrino SM 
gauge singlet held. 

Before proceeding, it is interesting to establish a parallel between the DS and the SM. In the SM, assuming only one 
generation and “turning off” the strong and the charged-current weak interactions, there are 7 massive Dirac fermions 
(m X 3, d X 3, e) and one massless left-handed Weyl fermion (the neutrino component of L, vl)- There is also a very 
large accidental classical global symmetry, 17(1)® (“green up-quark number,” “electron number,” “neutrino number,” 
“red down-quark number, etc). The charged-current weak interactions explicitly break this down to 17(1)® x U{1)(, 
for quarks (one for each color) and leptons, respectively. 

In summary, the chiral 17(l)i, model described here accommodates new stable massive fermions whose masses are 
proportional to a new mass scale, which we will associate with the dark matter. As we will discuss in Sec. El the 
new gauge interaction, combined with the interactions between the DS and the SM discussed in the next section, is 
sufficient to predict the dark matter relic density. In order to render the neutrinos massive, however, one is required 
to add new mediator fields charged under both the SM and U(l)y, and a new mass scale, described in ^Med- We deal 
with this issue in Sec. US 


III. THE MEDIATOR SECTOR, U(l)^ INTERACTIONS AND KINETIC MIXING 

We postulate the existence of heavy nonchiral degrees of freedom charged under the gauge symmetries of both the 
SM and the DS. For concreteness, let us choose a very simple possibility: one vector-like Dirac fermion X with mass 
A. X is a weak doublet with hypercharge -1-1/2 and 17(1)[/-charge —1 whose renormalizable Lagrangian reads (we 
define = (a;,®'^), where x,a:“ are left-handed Weyl fermions): 

Jfued = iXJ/>X + AXX - (^klLxc/) + Kaixff + h.c.) , (12) 


where 


= + (13) 

T“ are the Pauli matrices, and L'^ = is the SM lepton doublet. We will discuss most of the consequences of 

this part of the Lagrangian in the next section. Before doing that, we need to address the new gauge interactions. 

The fact that .if is invariant under 17(l)v x 17(1)[/ implies that one must also take kinetic mixing between the two 
17(1) field strengths into account: 


C/^ 7 TDliV p 

.^Kin-Mix —- 7 ^ 




(14) 


where represents the hypercharge field strength. In the next sections we will argue that phenomenological 
considerations impose a bound (/i/Sin ?7 < 10“®. We now show that this value is natural in the sense that quantum 
corrections do not destabilize it. 

If the kinetic mixing is set to zero at some scale, quantum corrections give rise to a nonvanishing 77 . These effects 
appear as a consequence of the mediator sector, which couples to both the SM and the DS 17(l)s, and generates a 
nonzero sin rj at the one-loop level. This one-loop correction to sin 77 is given by 


A sin 77 = C, 


9 9u 
IGtt^ 


log 


A 


(15) 




5 


where g' is the hypercharge gauge coupling, g is the renormalization scale and (7^ is an 0(1) coefficient that depends 
on the renormalization conditions. Taking for the combination C^g'g^, ^ 0.1, Eq. (1151) implies that sinry ^ 10“^ is 
natural. 

After electroweak and U{l)i, symmetry breaking, the physical gauge bosons A^, and are related to W^, 
Bfi by the following nonunitary linear transformation [l^ : 

( Bfj_ \ / cw — cos^svK ~ tanrysin^ sin^SvK “ tan?ycos^ \ / A^\ 

^/? = svy cos^cu/ -sin^cvy 

B^ J y 0 sin ^ sec ry cos ^ sec ry J \ J 

where cw and sw are, respectively, the sine and cosine of the weak mixing angle, while the angle ^ can be succintly 
written in the limit jy —>■ 0, Mz <C as (l^ 


sin^ ~-^svvsinTy. (17) 

Unless otherwise noted, we will mostly be interested in the limit in which Eq. dm holds. Notice that in this limit 
C < 7- 

In general, we find that DS particles acquire an 0{^) small coupling to the physical Z boson while SM particles 
acquire an 0{g) coupling to the physical Z gauge boson. These couplings are important because they represent the 
main SM-DS interaction and can be potentially probed experimentally. On the other hand, no coupling of the DS to 
the photon appears as it is expected since electromagnetism remains unbroken. Explicitly, the couplings of the DS 
particles to the Z are 


gv sin ^ - 2NlJ^Nl - ‘^XrI^XR + XL Qz^I^Xl) Z^ . (18) 


The SM particles also interact with the physical Z. In the same approximation as before - ry <C 1, Mz - 

the interaction can be written as 


e sin^Z^ 


tan 9w \ ^ + 


Ml 




Mi 


jif) 

cot 9w 


irf ^ - 


e sinry 
Cw 




(19) 


where the sum runs over all chiral SM fermions and and are the weak isospin and hypercharge of the fermion 

field / = UL,UR,dL,dR, _ For <C M| the term proportional to Y^^'> dominates. 

There is also mixing in the scalar sector. Again, assuming that the coefficient X^h of the X^h\(I>\'^\H\'^ term vanishes 
at tree level, a coupling X^h would be generated at the one-loop level of order 








IGtt^ 


log 



( 20 ) 


where we ignore a. Z — B one-loop diagram proportional to sin^ ^ (which can be considered a three-loop effect). We 
take the value in Eq. (EOl) to be the “natural” value for X^h- For small kr and Hi, X^h is naively more suppressed 
than sin?y. Moreover, as we will briefly argue in the concluding section, if one were interested in avoiding large finite 
corrections from the mediator scale A to the scalar masses-squared, one would be forced to impose KL,Ki < 10““^. 
Henceforth, we will assume for simpliciy that kinetic mixing effects are more significant than those related to scalar 
mixing. 


IV. NEUTRINO MASSES 

The mediator Lagrangian in Eq. m explicitly breaks the C/(l)jvi x U{ 1 )n 2 x U{1)^c x U{l)e down to U{1)l - 
what is normally referred to as lepton-number. That is, the introduction of new, heavy degrees of freedom necessarily 
renders all neutrinos and DS particles massive, but does it in a way that lepton number remains a good quantum 
number Q . Examples of recent models that make use of this idea are 0. 
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Upon integrating out the X field^, one generates the following dimension-five effective operator, assuming only one 
generation of SM leptons: 


-^(lVL)(i/<^) + h.c.. (21) 

It is important to note that the mediator sector we have chosen is an example. Any UV-completion that implements 
the symmetry breaking pattern described above would also lead to the operator in Eq. (EH- After symmetry breaking, 
this operator manifests itself as a Dirac mass between the 1+ fields and z/l, 


U KLKiVVA 

m, = 


A 


( 22 ) 


where v is the vev of the Higgs field, v ^ 10^ GeV. In more detail, after 1/(1)^ symmetry breaking, the three pairs of 
Weyl fields (z/l, 2j ) and ’ combine into three massive Dirac fermions, with mass matrix 

(Afi/)y = fijVtj,, (23) 


where /b = /y for i = 0,1, 2 and j = 1,2, while //q = mf for i = 0,1, 2. In the basis where dehned in Eq. ([T]), 

is diagonal, 

/ 0 0 \ 

M, = mf Ml 0 , (24) 

y 0 M 2 / 

where Mi, M 2 are the eigenvalues of naively of order v^. The masses m^, on the other hand, are parametrically 
smaller than Mi ,2 by a factor v/A. In the limit Mi ,2 ^ m^i 2 , the eigenvalues of AI^, are simply ttiq , Mi, M 2 . With 
respect to the SM weak interactions, the three massive Dirac fermions are a mostly active neutrino v and two mostly 
sterile neutrinos A^i, 2 . 

The fact that there are three generations of weakly interacting neutrinos in the SM, two of which are known to 
be massive, implies that the number of DS “families” Nf is at least two. Eor Nf = 2, there are three mostly active 
neutrinos, one of which is massless, and four mostly sterile states. All massive states are Dirac fermions. If all mostly 
active neutrinos are massive, Nf > 3 and there are at least six mostly sterile massive states. Eor Nf > 3, there is a 
mismatch between the number of vl and N states, and Nf — 3 antineutrinos remain massless. 

Experimentally, neutrino masses are known to be of order 10“^ eV or less which implies 


KLKpVtj, ^ 

A 


(25) 


The smallness of the neutrino masses can be attributed to A ^ or klKq I. Unless otherwise noted, we will 
assume that ^ v ^ 10^ GeV so if kqKl ^ 1, A ^ 10^^ GeV is required in order to “explain” the small neutrino 
masses. On the other hand, if A ^ 1 TeV, small neutrino masses can be “explained” if kqKl ^ 10“^°. This is, 
qualitatively speaking, not different from the standard Type-I seesaw mechanism Q. It is amusing to note that, in 
the limit kl = np and v = v^, the expression for small Dirac neutrino masses we obtain here, Eq. (I22L is identical to 
the one for small Majorana neutrino masses from the Type-I seesaw if we identify A with the right-handed neutrino 
masses and k with the neutrino Yukawa coupling. 


t We will mostly be interested in the physics at energy scales well below the mass of these mediator fields and could also introduce their 
effects by adding higher-dimensional operators to J(sm +.^03- We find that the introduction of a concrete model renders the discussion 
more transparent. 
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It is easy to diagonalize in the limit Mi, M 2 » mf*, mf*. M" = where 



( 1 

mf* /Ml 

m§/M2\ ^ 

Ul = 

-mf/Mi 

I 

0 +^(w) 


\ -mf/M 2 

0 

1 1 V W 



(26) 

(27) 

(28) 


In the case of we also only keep the leading order mf’/Mk terms. In the mass basis, the gauge interactions 

will couple the heavy, sterile states N to the active ly state. In more detail, the following interactions appear after 
diagonalization of 


rrii _ 

Wi^^ 


2cw 


-2g. 


m 2 _ 


2cu/ 


$ ~ ] N 2 L + h.c., 


(29) 


describing the potential decays Nk —>■ Zv, Bv. Moreover, from the standard weak couplings to the W^-boson, we find 


9 

V2 


■.D 




M^ ^ ^ 


N 2 L + h.c., 


(30) 


which allow for the interesting decays of the DS N particles into charged leptons, Nk —>■ . The strength of the 

coupling of A^s to the SM gauge bosons, proportional to the ratio between the neutrino masses and the masses of the 
sterile neutrinos, is a generic feature of these models. 

In the case in which the N particles are heavier than the weak bosons, the decay rate of, say, Ni —>■ iW is given by 




Gp[m?YMi 

87rV2 \MiJ 


(31) 


from which we obtain T ^ 0.4 s“^ for mf = 0.1 eV, Mi = 1 TeV. On the other hand, if the SM gauge 
bosons are heavier than the Ni fields, these decay to SM fields via off-shell H^-bosons and Z-bosons. In this case, the 
decay width of IV —>■ SM scales like Mf - a factor of Mf from kinematics times the “mixing parameter” squared, 
proportional to (mf /MiY'. For very light sterile neutrino masses - masses below 1 MeV - at the tree-level, only the 
N —>■ vvv decays are kinematically available, and the one-loop suppressed decay iV —>■ 1/7 also becomes relevant. The 
couplings associated to these decays are again proportional to /M. In summary, assuming all are of order 

the active neutrino masses, the sterile neutrino lifetimes range between tenths of milliseconds for Mjv 1 10 TeV to 
order 10® years for Mjy ^ 100 MeV and much longer than the age of the universe for lighter sterile neutrinos. 

It is intriguing that for sterile neutrino masses of order 10 keV, the active-sterile mixing angle squared (m^/M)^ 
is of order 10“^'^, in agreement with the recent “sterile-neutrino-as-dark-matter” interpretations of the 3.5 keV line 
Here, however, one needs to revisit the issue with some care since, in the early universe, the sterile neutrinos are 
kept in thermal equilibrium with the photons via flavor-diagonal Z and Z interactions, discussed in some detail in the 
previous section. These interactions determine their relic abundance, as opposed to the standard Dodelson-Widrow 
mechanism d, where active-sterile mixing determines the relic abundance of the mostly sterile states (for other 
possibilities see, for example, 0 ). We discuss early-universe related issues in more detail in Sec. [V] 

Sterile neutrinos can be produced in the laboratory mostly via their coupling to the Z and the Z [Z ^*'>, Z^*'^ —>■ NiNi), 
or the scalar fields associated to spontaneous symmetry breaking. The associated phenomenology and potential current 
constraints will depend on the mass of the sterile neutrinos but, for most masses, as discussed above, the W particles 
are effectively stable when compared to the time-scales of laboratory experiments and will manifest themselves as 
missing energy. 


§ We restrict the discussion to one generation of DS and SM fields, 
forward. 


The extension to three SM families and Nf DS families is straight- 








V. EARLY UNIVERSE COSMOLOGY AND DARK MATTER 


Here we provide a brief discussion of the thermal history of the universe described by the Lagrangian in Eq. O- 
As discussed in Sec. im in the scenario under investigation there are a number of accidental global symmetries that 
ensure the presence of stable particles. If there were no mediator states X, for every generation of new fermions, 
the corresponding fVi, N 2 would be stable, along with the massless “left-handed antineutrinos” v'^. Furthermore, yi, 
the lightest of the Xi fields, is also stable. The mediator interactions discussed above, other than rendering the Ai 2 
unstable on cosmological time-scales do not play a significant role at temperatures much smaller than A. In what 
follows, we estimate the current constraints on the DS from cosmological measurements, and discuss whether the y 
particles can explain the dark matter puzzle. For concreteness, we will use the following values for the DS parameters 
in our estimates: 


Ml - M 2 - 1 TeV, M^ - 5 TeV, M^ - 500 GeV, 5 ,. - 0.1. (32) 

We denote the temperatures of the SM and DS plasmas by T and T respectively. The distinction is only relevant 
below a certain temperature Tdec at which the interactions that couple the two sectors fall out of thermal equilibrium. 
Above Tdec there is a single plasma composed of SM and DS states. Also, let g* represent the effective number of 
relativistic degrees of freedom at temperatures greater than Tdec while g* and g* represent the degrees of freedom 
below Tdec for the SM and DS plasmas respectively. To simplify things we assume that no reheating takes place at 
the time of decoupling so that at Tdec we have the matching condition: 

gdec=5dec+53ec- (33) 

The thermal history of this model proceeds as follows [llj . As the universe cools down, DM becomes nonrelativistic 
and freezes out eventually, leaving a thermal relic of yi particles. Annihilations of the heavy, stable yi particles 
into lighter or Ni dominate over annihilations into SM states and are the main processes that determine the relic 
density. Shortly after freeze-out, all the states in the DS are nonrelativistic except for the These remain coupled 
to the SM plasma up to Tdec- From this point on, T is no longer necessarily equal to T and the two gases must be 
treated independently. After the time of decoupling, adiabaticity of the evolution of the two sectors then imposes the 
constraint 


ffSec 

9*T^ 9Lc ■ 


(34) 


We start by estimating Tdec and how it is constrained by observations. Given the choice of parameters Eq. dsa), this 
occurs at temperatures well below the electroweak and U{l)i, phase transitions. At temperatures T = T < 100 GeV, 
the only relativistic state in the DS plasma are the The presence of new, light particles in thermal equilibrium 
with the SM gas is contrained by a variety of observations. This is true of the antineutrino fields which are 
effectively massless. These degrees of freedom persist throughout the thermal history of the universe and contribute 
to its expansion rate, especially when the universe is radiation dominated. More concretely, they contribute to the 
parameter Nes, which parameterizes the number of relativistic degrees of freedom: 


Veff = + AVeff , AVeff = ^ . (35) 

Pi/ 

Here is the SM contribution while ANgg comes from new physics and is defined as the ratio between the energy 
density of extra radiation peR and that of one SM neutrino p,^. 

The SM particle content and interactions translate into = 3.045. Measurements constrain iVefi at the time 
of big-bang nucleosynthesis (BBN) and at the surface of last scattering. In particular, if the gas of states were at 
the same temperature as the active neutrino gas at the time of BBN, they would contribute ANgg = Ng Since we 
need at least two DS generations in order to account for two massive neutrinos, we have Af > 2, a possibility that 
would be in tension with the current bounds if it were to translate directly into ANgg. Measurements of AN^g come 


^ For masses of the sterile neutrinos Mi, M 2 ~ 1 TeV, their lifetimes are a fraction of a second and hence “safe” with respect to 
constraints from big-bang nucleosynthesis, since the injected particles have sufficient time to thermalize. We assume therefore that the 
sterile neutrinos are mostly harmless as far as cosmological observables are concerned and ignore them henceforth, unless otherwise 
noted. 
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from BBN and the cosmic microwave background (CMB). The Planck collaboration has recently published values 
Ngff = 3.15 ± 0.46 at 95% C.L. [2^ while in [2^ the bound < 1.5 at the time of BBN was found. 

On the other hand, if the DS gas decouples early enough from the SM its temperature T at the time of big bang 
nucleosynthesis can be significantly lower than that of the active neutrinos T. If that is the case, 


A A. ^9* 

AiVeff = 


4r 

7 V^dec^y 


(36) 


where we have used Eq. dal. Assuming Nj = 3 generations in the dark sector g* = = 6, and taking for g* 

its value at the time of neutrino decoupling g* = 10.75, we find that in order to satisfy AA^eff < 1, it is enough 
to have g^^^ > 27. Hence, it suffices that the decouple before the QCD phase transition, which corresponds to 
9dec ~ 61.75, in order to satisfy the bounds comfortably. In that case we find ATVeff = 0.33, in agreement with Plack 
bounds. 

For later use, it proves useful to define the ratio between the DS and SM temperatures today (indicated by the 
subscript ‘0’) which has remained constant after the photon reheating by electron-positron annihilation:. Using 
Eq. dSl we find: 


T 

ro = ^< 0.56. (37) 

To 

Imposing that the DS should decouple from the SM before the QCD phase transition translates into bounds on 
the couplings of the and SM fermionic currents to the Z and Z bosons respectively. We note that thermal 
equilibrium with the SM is very similar to the thermal equilibrium of active neutrinos since these also interact with 
the rest of the SM gas via weak interactions. Since the decoupling temperature of active neutrinos is proportional to 
, we roughly estimate the decoupling temperature of the from the SM gas, for Mz M^, as follows: 


rdec ^ J 

where Gv has been obtained from Eqs. ()18I19L 


/ Mzgv 



G 


sin^ C gl 

Ml 


(38) 


(39) 


Requiring that the decouple from the SM gas before the QCD phase transition, j'dec^j'dec < q .^g obtain 
sin?? ^ 10“^. This matches the discussion in Sec. IIIIl where we argued that sin?? ^ 10“^ is natural (see Eq. (fTSl) '). 

We proceed to estimate the DM relic density given the values of the parameters listed in Eq. (ESI). Define as usual 
the time variable x = and let Xfo be its value at the freeze-out temperature (we will show that Tfo = Tfo)- Xfo 

can be obtained from the Boltzmann equation and is typically of order 

Xfo (40) 

^ fo 

It is well known that the contribution of xi to the energy budget of the universe today is inversely proportional 
to its thermally-averaged annihilation cross-section, (cTann??)- For s-channel annihilations we have: 


(CTann??) ^ C^iXl + Xl ^ ^ NfQ' 


2 9Z 




(41) 


where = 0(10) was defined in Eq. dill) and Nf is the number of dark flavors. Putting in some numbers we find 


(o-ann??) -- 3 X 10 pb X 


1 TeV 
M. 


XI 


NfQ 


20 


(S'- 


(42) 


For a DM number density , the DM fraction today is defined, as usual, as 


= 


PxO _ f OooUoo\ _ My,^Yoof^ ( ttooToo \ (OdecTde 


00^0 


Per y UdecTdec j y ^0^0 


MxrYoof^ ( aooTA 


y ??dec^d e 




Pc 


s; 


( 43 ) 
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where Y = /T^ and the c» subscript refers to its asymptotic value. In the above equations we have used that 

Tdec = 2dec and T^o = Too- Also, notice that no reheating event occurs in the DS after T^ec and hence OdecTdec = clqTq. 
A good approximation for Yoo is 


XioH{x = 1 ) 

M^^{aannV) 


(44) 


where the Hubble rate at a; = 1 is given by 

= 1 ) = ° 

Plugging this into Eq. (l43l) we obtain 


4n^GQ*{x = 1 ) 

45 


Mi 


^xo 


47r3G'g*(x = 1) 


45 


(o'ann'c)Pcr 9^ 


(45) 


(46) 


Numerically, 




10-2 

annn) 


pb. 


(47) 


Given the value for the thermally-averaged cross-section, estimated above (Eq. (H^ l. can be made to agree with 
the cold dark matter contribution to the energy budget of the universe, = 0.1188 ± 0.0010 [^ . 

The xi particles scatter off of ordinary matter via both Z and Z exchange, both couplings (i.e., the coupling of 
Xi to the Z-boson and that of the SM fermions to the Z-boson) suppressed by the small kinetic mixing parameter. 
The cross-section for xi~nucleus scattering is, in the limit that the dark matter is much heavier than the scattered 
nucleus and M| Mf [^** 


ct(Xi + N ^ Xi + N) 


^ mlf ggsin^^Q^ 

sin^ 9w cos2 9w 


(1 — sin^ 9w)^Z^, 


(48) 


where a is the fine-structure constant, Qy = {4+ [(3 -I- 2|V;,^5|2])2, 49 < Qy < 64 is the square of two times the U{l)i, 
vector-charge of the dark matter candidate and mjv is the mass of the nucleus with atomic number Z. It is interesting 
to note that in the usual WIMP scenario, the direct detection cross section is dominated by scattering off neutrons 
while in this case, the proton contribution is the most relevant one. For xenon, the cross-section per nucleon, defined 
as cr^p = (j{xi + Xe —>■ Xi -I- 'Z.e)rri^/m\^A^^ where rup is the nucleon mass, is 

a^p = 1.4 X sin2 ^ ^ ^ 

For = 5 TeV, the LUX experiment constrains a^^p < 6 x 10“^^ cm^ at the 90% confidence level [2^, which 
translates into 

“xl»-'x(^)(^). (50) 


The next round of direct-detection experiments, assuming that WIMPs are not detected, will start to seriously 
constrain xi the dark matter. The estimates for the relic density and constraints from direct detection are in 
agreement with more general results for “electroweakly coupled” DM [Z7l| . 

Finally, the model also predicts signals for indirect detection experiments. In regions where the density of xi 
particles is large, they can annihilate into light SM particles, right-handed neutrinos, or sterile neutrinos. In the limit 
discussed above. My, ^ Mz,M^, we expect the and, if kinematically accessible, NN final states to dominate. 
Right-handed neutrinos are virtually invisible. At high energies, the states are, for all practical purposes, massless, 
and interact with ordinary matter via Z and Z exchange, both cross-sections suppressed by sin^ ^ relative to that of 
ordinary high energy neutrino-matter scattering. The NN final states are more interesting, given that the N particles 
will decay into high energy neutrinos and charged-leptons, as discussed in Sec.H^ Hence, the process xiXi is 

expected to yield high energy (energies < Myi) electrons, positrons, and neutrinos. 


Even if one is interested in the limit M? ^ .Z-boson exchange is still comparable to Z-boson exchange in the limit where the 

momentum transfers are much less than as discussed in Sec. nni 
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VI. ABELIAN GAUGE THEORIES WITH CHIRAL FERMIONS 

In this section we describe a general method for producing models of a Dark Sector with a fermion content that is 
chiral and anomaly-free (i.e. models in which there does not exist a right-handed partner with opposite charge for 
every left-handed particle). This issue was also studied in detail in [T^ . Some of the results presented here were also 
discussed in [l^ . 

Let be a simple Lie group and H be its Cartan subalgebra. Any H gH defines a 17(1) subgroup of The states 
in an arbitrary representation r transform under this 17(1) with charges given by the eigenvalues of 17^, the Cartan 
generator written in the representation r. 

A sufficient condition for a gauge 17(1) to be anomaly free is to impose that charges and field content are in 
correspondence with the eigenvalues of Hr for a group 1# that is anomaly-free. It is well known that the classical 
groups SU{2), SO{n) for n ^ 6, Sp{2n), G 2 , F 4 , Ej and E% are all anomaly-free. For ^ any of these groups, a choice 
of representation r and Cartan generator H determines a 17(1) gauge theory that is anomaly free. 

If we demand such gauge theory to be also chiral, the set of possible groups and representations is reduced consid¬ 
erably. All of the above groups except for those of the form SO{4:n + 2) only have real representations. For a real 
representation, the matrices Hr and —H* are similar, Hr = U{—H*)U'^ for some U in 1#. Hence, the eigenvalues of 
Hr and —H* must be equal. Since Hr is hermitian, its eigenvalues are real. From here we conclude that if A is an 
eigenvalue of Hr, then —A also is. Hence, real representations do not lead to chiral models. 

In order to generate anomaly free chiral models one needs a group with complex representations that is anomaly 
free. «S'O(10) is the smallest such group since S'0(6) is anomalous. The smallest complex representation of «S'O(10) is 
the spinorial 16. A generator belonging to the S'O(IO) Cartan subalgebra in this representation has the general form 

H (a, b, c, d, e) = —diagja + h-\-c + d + e, —a + b + c + d — e, a — b + c + d — e, —a — b + c + d + e, 

a + b—c + d— e, —a + b — c + d + e, a — b — c + d + e, —a — b — c + d — e, 

a + b + c—d — e, —a + b + c—d + e, a — b + c—d + e, —a — b + c — d — e, 

a + b—c—d + e, —a + b — c—d — e, a — b — c—d — e, —a — b — c — d + e} (51) 

where a, b, c, d and e are arbitrary real numbers and V is a normalization factor. Two cases can be distinguished: 

• At least one of a, b, c, d or e vanishes. Then, —H also belongs to the Cartan subalgebra. 

• None of a, b, c, d or e vanishes. In this case, —H is not a generator in this representation. 

In any case, since iS'O(lO) is anomaly free, we have 

Tr [i7] = Tr [H^] = 0 , (52) 

as promised. That is, as long as none of a, b, c, d or e vanishes, the elements in the diagonal above are chiral solutions 
to the anomaly equations 


16 16 

i=l i=l 

where the qi are U{1) charges with respect to a 17(1) gauge group. 

Without too much loss of generality, we restrict the rest of our discussion to integer charges. From this requirement, 
either none, two, or four out of a, b, c, d and e can be half-integers, the rest must be integers. Now, notice that 

H{—a, b, c, d,e) = H (a, —b, c, d,e) = H (a, b, —c,d,e) = H (a, b, c, —d, e) = H (a, b, c, d, —e) = —H{a, b, c, d, e). (54) 

Hence, up to an overall minus sign in all charges, it is enough to consider a, b, c, d, e strictly positive. The largest 
U{1) charge in this case is given by a + b + c + d + e and notice that it appears only once. The absolute values of 
all other charges is necessarily smaller. Hence, for any anomaly-free model derived in this way, there can be only one 
state with the highest 17(1) charge. In particular, this state cannot be charged under any extra nonabelian gauge 
symmetries that commute with the 17(1) in question. In the SM, that state is the left-handed anti-electron e'^, which 
has hypercharge 6 in units of the smallest known hypercharge, that of the doublet quark field Q. Note that e'^ is not 
charged under color SU{3) or electroweak 517(2). 

The following values 


a=b=c=d=- 


e = 1, 


(55) 
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yield a solution for the anomaly equations with the smallest highest charge possible, equal to 3 units in the normal¬ 
ization defined above, in which all charges are integers. Discarding nonchiral pairs (i.e., charges which are equal in 
magnitude and opposite in sign) and vanishing charges, this solution has the following particle content: 

3x1, -2x4, 1x5, (56) 

where n x to means to particles of charge n. All together there are ten chiral fields. It can be quickly checked that 
the anomaly conditions are satisfied. 

With highest charge equal 4 we can take the a, 6,c, d, e combinations {1/2,1/2,1,1,1}, {1/2,1/2,1/2,1/2, 2} or 
{1/2,1/2,1/2, 3/2,1}. The last of these yields a particularly simple particle content (only 7 states): 

4x1, -3x3, 2x2, 1x1. (57) 

With highest charge equal to 5 an even simpler particle content is obtained (6 fields): 

5x1, -4x2, 1x3. (58) 

Then, what is the minimum particle content possible that is chiral and allows for at least one anomaly-free 17(1) 
gauge symmetry? The following solution of Eqs. (I53|) with highest charge 10 requires only five fields: 

10x1, -9x1, -7x1, 4x1, 2x1. (59) 

On the other hand, there are no solutions to Eqs. (1551) with three fields due to Fermat’s Last Theorem. It is also 
relatively easy to prove that for four fields only nonchiral solutions exist^b Hence hve chiral fields is minimal. 

The model we explore in the previous sections is not minimal - ten fields, highest charge 5 - according to the two 
criteria defined above. It, however, still “more minimal” than the SM, which has hfteen chiral fields per generation 
and highest charge 6. 

The model associated with Eq. (l56l) . if one assumes a Higgs-like field with charge -1-1 and a nonzero vacuum 
expectation value, would, in general, contain four Dirac fermions and two massless chiral fermions (linear combinations 
of the one 3 and five 1 states). Similarly, the model associated with Eq. (l57l) . if one assumes a Higgs-like field with 
charge -1-1 and a nonzero vacuum expectation value would, in general, describe three massive Dirac fermions and a 
massless chiral fermion, 1. This model is similar to the one described in detail in the previous sections, minus the 
Dirac “sterile” neutrino states. 

It is amusing that in the case of the smallest model (five fields), associated to Eq. (1551) one needs more than one 
Higgs-like field in order to render more than one pair of chiral fermions massive. If one were to add scalar fields with 
charge -1-1 and -1-6 one would end up with two massive Dirac fermions and one massless state, —7. 


VII. COMMENTS AND OUTLOOK 

The fermion content of the SM is chiral and before EWSB, all masses are identically zero. Only after the spontaneous 
breaking of the gauge symmetry does the theory allow for massive fermions and all masses are proportional to the 
same mass-scale. A side-effect of the chiral nature of the SM fermion content is that some fermions end up massless 
(i.e., the neutrinos) and some of the massive objects in the theory are, thanks to accidental global symmetries, stable 
(i.e., the proton). Inspired by the SM, we posit that all “light” fermions must be chiral as far as some gauge symmetry 
is concerned and therefore massless for as long as the symmetry is manifest. Masses for all such fermions should come 
as the result of a Higgs mechanism that spontaneously breaks the gauge symmetry at some energy scale. 

The scenario described above invites two main questions. First, can one construct models within this paradigm 
generically and systematically? How? Second, are these models useful to address outstanding puzzles in fundamental 
physics, including the dark matter puzzle and the origin of nonzero neutrino masses? In the bulk of this paper, we 


tt We write Eq. 1531 1 in such a way that all qi are positive. Then there are essentially two cases: 


1. qi = q 2 -\-qz-\- q^. Hence 

ql = (92 -I- 93 -I- 94 )® >92+93+94 

so there can be no solutions in this case. 

2. + 92 = 93 + 94 Q-nd 5 ^ + 91 = 91+ q^. Then, one must also have 9192 = 9394 - Any solution to this set of equations involves two 
pairs of equal charges. 
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dealt with the second question. We fleshed out a model of a dark sector (DS), consisting of a new, spontaneously 
broken along with chiral fermions charged under U{1)^. This model indeed leads to a candidate DM particle in 

the form of a new WIMP-like massive stable state. A relic abundance for that state appears as a result of a standard 
thermal freeze-out mechanism occuring in the DS. We showed that the breaking scale can be chosen such that 

the relic density of the WIMP-like object is consistent with the evidence for dark matter in the universe. 

The sample-model also contains massless chiral fermions, which as far as the SM gauge group is concerned, look 
like right-handed neutrinos. Nonzero neutrino masses require the introduction of a new, potentially very heavy, 
mediator mass scale that connects these states with the massless left-handed neutrinos of the SM. We discussed in 
detail an ultraviolet complete scenario that preserves lepton number but nonetheless leads to parametrically small 
Dirac neutrino masses. anomaly constraints also require the existence of at least four mostly sterile Dirac 

neutrinos - two sterile neutrinos are required per generation of new fermions, and the neutrino data require at least 
two of these generations - that mix with the SM neutrinos. The mixing between active and sterile neutrinos is related, 
qualitatively, to the observed nonzero active neutrino masses. 

Some extra features of the model might be considered appealing. The scalar sector of the DS is minimal in the 
sense that one field is enough to give mass to all the DS particles except for the RH neutrino. Moreover, the fermion 
content is also minimal in the sense that only two states are stable on cosmological scales, the DM and the RH 
neutrino. Thus, the only stable particle that the model adds to the SM ones is the DM fermion. All in all, this model 
is rich enough to provide DM, neutrino masses and some interesting phenomenology, and yet it is simple enough not 
to modify the SM dramatically. 

For the sake of concreteness, we brushed aside several possibilities that may be worthy of further exploration. For 
this model, we chose a mediator sector that consisted of a new vector-like fermion charged under both the SM and the 
C/(l)jy gauge symmetries. A different mediator sector would lead to other sources for nonzero masses for both the SM 
and the DS “neutrinos,” including some that would render all of these particles massive Majorana fermions. Different 
possibilities can be explored qualitatively by appreciating that there are two dimension-five operators other than 
Eq. (I^Tl) one can construct, given the particle content and the gauge symmetries discussed in detail here. Explicitly, 
one may have {LH){LH), the so-called Weinberg operator, and (l+0)(l+(/)). The former leads to nonzero Majorana 
masses for the active neutrinos, proportional to /M^ while the latter leads to “right-handed neutrino” masses 
proportional to where is the generalized lepton-number-breaking scale. This possibility was considered 

recently in The scale is, in general, not related to A, the mass of the mediator field X. If ^ A, the results 
discussed here remain valid. This would happen, for example, for ^ 10^® GeV, the grand-unification scale. In 
this case, both left-handed and right-handed neutrinos receive nonzero Majorana masses, but these are much smaller 
than the Dirac masses and the neutrinos would be pseudoDirac fermions that, for most experimental purposes, “look 
like” Dirac neutrinos. 

More speculatively, while we extended the SM gauge symmetry with a new [7(1 )i,, more complicated DS gauge 
symmetries are also a possibility worthy of pursuit. In particular, it is easy to check that in the model discussed in 
detail in this paper, one can add a new, nonanomalous^^ SU{2)y gauge symmetry, under which 1^), {2^, 2^), 
(3(_,3?.), (4j,_,4^) all transform as doublets. This SU{2)^ is not spontaneously broken by since </> is an SU{2)y 
singlet. There would be, nonetheless, four massive Dirac fermions per DS generation, pairwise degenerate due to 
the unbroken SU(2)i, symmetry. At low energies and for a small enough number of dark matter generations, the 
SU{2)i/ interaction would confine and those states charged under it would manifest themselves as dark “mesons” or 
“baryons”. In particular, the lightest, stable such state could be the dark matter§§. Note that this SU{2)i, cannot be 
naively identified as the SM SU (2 )l since, if this were the case, all members of DS doublets would acquire electric 
charge ±1/2 after EWSB. 

It is interesting to enquire whether the model proposed here is natural, at least as far as it is defined in, for example, 
[^[^. In more detail, naturalness in this context translates into requiring that finite loop-corrections to the masses 
of H and do not supersede the masses themselves. The presence of the messenger field A, for example, can lead to 
very large corrections to the masses of both scalars. At the one loop-level, naturalness translates into, roughly, 

^A^, ^A^ < (100 GeV)^ (60) 

OTT^ OTT^ 

According to this criterium, keeping in mind that the active neutrino masses, of order 10“^ eV, are proportional to 
the theory is natural - in the sense that the A mass scale does not destabilize the electroweak and U{l)y 


As far as the SU{2)iy is concerned, the fermions are vector-like. 

The case where the SU{2)v is infrared free is probably ruled out by the fact that long-range dark matter self-interactions are excluded 
by observations. 
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breaking scales - as long as A < 10^ GeV (assuming ~ kq and v ^ ^ 10^ GeV). 

Independent from the mediator sector, the presence of two mass-scales, v and is “unnatural” [H, [sOl unless 

V ^ or if the coupling X^h between the two sector is very small. Naturalness from the neutrino sector naively 
requires, as discussed above, < 10“"^, and one can hence argue that tiny X^h values are, in some sense, expected 
(see discussion in Sec. uni and Eq. (1201) 1. This implies that even if one takes naturalness considerations into account, 

V can be significantly smaller (or larger) than 

Many other potential consequences of the model were not discussed here. Since the neutrinos are Dirac fermions, 
standard leptogenesis does not work. It is, therefore, reasonable to ask whether the scenario discussed here can 
accommodate a baryogenesis mechanism. For example, would Dirac leptogenesis work? We also did not consider 
the possibility that the U{1)^ gauge boson is light (M^ < 1 GeV) or very light {M^ ^ 1 MeV), nor did we explore 
the consequences of postulating that the dark matter is light. We also did not consider the possibility that the sterile 
neutrino states - whether they be hot, cold, or warm - could make up most or part of the dark matter. 

Finally, we return to the first question posed above: how to construct such models generically and systematically? 
In Sec. EH we discussed a simple yet powerful way of generating an infinite number of chiral, anomaly-free U{1) gauge 
theories, of which our model is only one example. Each model found by this method will have its own idiosyncrasies 
and many should allow one to address the current puzzles in particle physics. It also seems more-or-less straightforward 
to use the same S'O(I0)-inspired mechanism in order to identify anomaly-free theories with chiral fermions that are 
invariant under larger gauge symmetries, like the SM. 
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